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ENERGY CONCENTRATION AND A PRIORI ESTIMATES FOR B2 
AND G2 TYPES OF TODA SYSTEMS 

CHANG-SHOU LIN AND LEI ZHANG 


Abstract. For Toda systems with Cartan matrix either 82 or G2, we prove that 
the local mass of blowup solutions at its blowup points converges to a finite set. 
Further more this finite set can be completely determined for 82 Toda systems, 
while for G 2 systems we need one additional assumption. As an application of 
the local mass classification we establish a priori estimates for corresponding 
Toda systems defined on Riemann surfaces. 


1. Introduction 

Let {M,g) be a Riemann surface with area equal to 1 and Gauss curvature equal 
to a constant Kq. We consider the following Liouville equation 

N 

(1.1) Agu + 2{e‘‘-Ko)=4KY,aj6pj in M 

7=1 

where is the Laplace-Beltrami operator (— A,> 0), ay > — 1 and 5pj is the Dirac 
measure at p, G M. The geometric meaning of (11.11) is that for any solution u, the 
new conformal metric ds^ = e"g has the constant Gauss curvature equal to 1 outside 
the singular points {p/}. Near each pj, using a complex coordinate z satisfying 
z{pj) = 0 we have = C?(l)|zp“^ for |z| near 0, hence , the new metric ds^ is 
degenerate at py and is called a metric with conic singularity. Equation (11.11) and its 
general form, the so-called mean field equation, have been extensively studied for 
many decades. For example see |I7][T0l[TT][37l[39lllQl|43l and the reference therein. 
In particular, for the case M being the standard sphere or a torus, the equation can 
be written as ( replacing u in (11.11) by « + log 2) 

( 1 . 2 ) AgU + e‘'= 47i'^aj5pj. 

j 

Recently equation (11.21) is found to have deep connection with the classical Lame 
equation and also the Painleve VI equation. For example for the Painleve VI equa¬ 
tion with certain parameters, some non-existence theorem of (11.21) plays a key part 
in the proof of the smoothness of the solutions with unitary monodromy group. 
The interested readers may read into f91 and |[T2l for more in-depth discussions. 
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The natural generalization of (11.11) is the so-called Toda system 

n n 

(1.3) ^gUi + £ = 471 ^ aij5p,. 

;=i ,/=i 


where kij is known as the Cartan group of some simple Lie algebra. For the case of 
A„, type (see explanations below), the coefficient matrix A = (a,/) is expressed by 


(1.4) 


/ 2 -1 0 ... 0 \ 

-1 2 -1 ... 0 



-1 2 -1 
0 - 12 / 


It is well known that the general Toda system (11.31) is closely related to geome¬ 
try 15] 13 [Ml 1221 and the gauge theory in many physics models. For example, to 
describe the physics of high critical temperature superconductivity, a model of rel¬ 
ative Chem-Simons model was proposed and this model can be reduced to anxn 
system with exponential nonlinearity if the gauge potential and the Higgs field are 
algebraically restricfed. Then the Toda system (11.31) with (11.41) is one of the limit¬ 
ing equations if the coupling constant tends to zero. For extensive discussions on 
the relationship between Toda system and its background in Physics we refer the 
readers to l4l[l6l[l7][Il]|2Tl|25l|26l|27l|35l|43l and the reference therein. 

If the rank of the simple Lie Algebra is 2, there are three types of corresponding 
Cartan matrices of rank 2: 


A2 = 


2 

-1 



B2{=C2) 


2 

-2 



G2 


2 -1 
-3 2 


These rank 2 matrices are the simplest examples of their more general forms. 
In general, there are four types of simple non-exceptional Lie Algebra: 
and Dm whose Cartan subalgebra are sl{m + 1), so{2m + 1), sp{m) and so{2m), 
respectively. Corresponding to each of the four types of Lie Algebra there is a 
Toda system. Solutions of Toda systems are closely related to holomorphic curves 
in projective spaces. In particular, from the classical Pliicker formula we see that 
any holomorphic curve gives rise to a solution of the Am type Toda system and the 
branch points of these curves correspond to the singularities of the solutions. On 
the other hand, if we integrate the A„j Toda system, a solution defines a holomor¬ 
phic curve in CP” af least locally. The interested readers may see ||291 for further 
discussions of this respect. 

In llT3l and PTl it was noticed that the equation (11.31 ). like (11.11 ). is also an 
integrable system. The integrability has been further discussed in |j29l . As we 
mentioned earlier, equation (11.31 ) has deep connections with algebraic geometry, 
modular forms and the Painleve VI equation. So it is natural for us to study (11.31 ) 
from analysis viewpoints as well as the perspectives from integrable systems. From 
the analytic side, the most important issue is to derive a degree counting formula 
for equation (11.31 ). a generalization of the previous works of Chen-Lin lITOl fTTl 
for (11.11 ) and Lin-Zhang 1(31] l32l l33l for general Liouville systems. However this 
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generalization is very challenging because the bubbling phenomena are more com¬ 
plicated and the concentration has not yet been proven even for SU (3) Toda system. 
So far in this direction the degree counting formula has only been proved for the 
simplest case of SU (3) Toda system, see Il28]l . 

In this paper we initiate the analytic program for the systems B 2 and G 2 . The 
main purpose is to establish the a priori bound for non-critical parameters. Similar 
to (11.11) we consider the Toda system of the mean field type: 


(1.5) 


AgUi 


n 


y=i 


hje^j 


1) = 0,i = I, ..,n 


where h[,...,hn are positive smooth functions on M, pi,...,p„ are positive con¬ 
stants, and the solutions are in the space 

//h” = {v = (vi,...,v„); f VidVg = 0,VieL\M),yvieL\M) i=l,..,n}. 

Jm 

For B 2 or G 2 , the coefficient matrix A is 


A = 


2 

-2 



or 



-1 

2 


, respectively. 


One of the most important and intriguing issues of Toda system in general is the 
blowup phenomenon. A point p is called a blowup point if, along a subsequence, 
a sequence of solutions {u^ = satisfies 

max max nf = maxMf(p^) ^ ooj p, 

i B{p,S) i 


Undersfanding fhe asymptotic behavior of blowup solufions near ifs blowup poinf 
is crucial for many imporfanf quesfions related to fhe Toda systems such as a priori 
esfimafes, degree counting formula, exisfence resulfs and mulfiplicify resulfs, efc. 
We say a sequence of blowup solutions {u^} possesses energy concenfrafion if for 
some i, max m* 00 , g"; tends to a Dirac measure as ^ 00 . The purpose of this 

article is to study the blowup phenomenon and the energy concentration for (11.5b 
with B 2 or G 2 matrix. For simplicity we assume 


2 

(1.6) Aui -|- ^ Uijhje'*^ = 0, in B(0,1), / = 1,2 

;=i 

where B(0,1) is the unit ball in (throughout the paper we use B{p,r) to denote 
the ball centered at p with radius r), A = (a;y) 2 x 2 is a B 2 or G 2 matrix. 

In more precise terms we let = (mj , be a sequence of solutions to 

2 ^ 

(1.7) Au'l + Y,aijh]e"i =0, in B(0,1), / = 1,2 

2=1 

where h \, are two sequences of positive smooth functions with uniform bound: 

(1.8) ||/rf||c2(fi(o,i))<C, inB(0,l), V/ = l,2, 
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the origin is the only blowup point for = {u\,u\), which has bounded energy 
and oscillation on 5B(0,1): 

' max,maxB(o,i) u\ max,maxjfc:cB(o,i)\{o} “f < C(K), 


(1.9) 


u\{x) — u\{y)\ <C, Vx,y G 5B(0, 1), /=1,2. 


I 4(0,1)/=1,2. 

The main purpose of this article is to study the following quantity: 

(1.10) a, = lim lim/=1,2. 

S^0k^«>27l Jbs 

Note that the inside limit: limj;^oo, which is to be taken first, is understood as 
taken along a subsequence of with the same notation. The second limit: lim^^o 
indicates that we consider how the energy (the integration of ) concentrates 
at 0. 

Note that the oscillation finiteness assumption in (11.91) is natural and generally 
satisfied in mosf applicafions. The energy bound in (11.91) is also nafural for sys- 
term/equafion defined in two dimensional spaces. 

The first main result is: 


Theorem 1.1. Let A be the B 2 matrix, the blowup solutions = (nf,n^) o/dIZD 
satisfy ( 17.9D . and ( 17.81) holds for h^ = (/jj,/i|). Then ( 01 , 02 ) defined by ( 17.701) is 
one of the following types: 

(2,0), (0,2), (4,2), (2,6), (4,8), (6,6), (6,8). 

As an application of Theoem 11.11 we consider the B 2 Toda system defined on 
a compacf Riemann surface {M,g), whose volume is assumed to be 1 for conve¬ 
nience. 

AgU\ +'^P\( ~ 4 “ J„h2e‘‘2dVg “ 1 ) = 0 ) 


( 1 . 11 ) 


hie"l 


■ - 1 ) + 202 ( 77 - 4 (^* 


. 4“2 ^Plij^hie-idV, 

where /ii ,/i 2 are positive smooth functions on M, Ag is the Laplace-Beltrami oper¬ 
ator {—Ag > 0), Pi,P 2 are positive constants. 

It is well known that for solutions of the following general system on M: 


AgUi -\- aij(^p 

7=1 


hje'^2 


^ fMhje^jdVg 


- 1 )= 0 , / = 


the corresponding variational form is 


^p(“) = T L - L Pi log /1 

i j—I 7—1 Af 


hje'^idVg. 


where {d^)nxn = { aij)nL - 

For solutions of (11.111) in we prove the following a priori estimate: 
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Theorem 1.2. Suppose pi,p 2 > 0 and none of them is equal to a multiple of 4k, 
then 

l^iWI / = 1,2 

where u = {ui,U 2 ) is a solution of Al.ll\) in 

Next we consider the locally defined G 2 Toda system: 


( 1 . 12 ) 


AMj' + 2/i^e“i=0, 


[ - 3h\e‘^' + 2/j|e“2 = 0, in B{0, 1) 

where the assumptions on u^ = {u\,U 2 ) and h^ = are the same as in the B 2 

case. Then in this case we have 


Theorem 1.3. Let ( 01 , 02 ) be defined for blowup solutions u^ of as Al.lOk 
Suppose di.9l) holds for u^, then ( 01 , 02 ) satisfies 

3 af + 0I - 30102 = 601 + 202. 

If in addition we have 

lim — [ < 4 + 2 V 2 , lim — [ h\e'*^ < 10 + 2\/7, 

Ik .Jb(o,i) 2k .Jb(o,i) 

(01,02) is one of the following types: 

(2,0),(0,2),(4,2),(2,8),(4,12). 

As an application of Theorem 11.31 we consider the following G 2 Toda system 
defined on the compact Riemann surface (M, g): 


(1.13) 


Ag“l +2pi(j^/,j‘g»lrfy^ 1) P2{ j^hlgU2dVg 1)“0> 




A.gU2 3 Pi ( 1 ) + 2 p 2 ( 

Jj^h2e'^2dVg 

Comesponding to Theorem 1 1.21 we have 


Theorem 1.4. Let u = (ui,U 2 ) be a solution of di.iJI) in i7pi€(0,(8 + 

4 s/2)k),P2 G (0, {20+4s/1)k) and none of them is equal to a multiple of 4k, there 
is a constant C independent of u such that 

|m/| < C, / = 1,2. 


Theorem 1 1.2 1 and Theorem [T4] make it possible to calculate the Leray-Schauder 
degree dp^p^ for (pi,p 2 ) in intervals determined by multiples of 4?! and the unrea¬ 
sonable energy upper bound in (11.131) . In order not to make this article exceedingly 
long we shall address the degree counting formula in a subsequent paper. 

To end the introduction we describe the outline of the proof of Theorem 11.11 
and Theorem 11.31 First we use a selection process to determine a finite number of 
mutually disjoint bubbling disks. The idea of the selection process was first intro¬ 
duced by Schoen |[3^ and is very useful for prescribing curvature type equations. 
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In if^ Lin-Wei-Zhang applied this method to locate huhhling disks for systems 
of equations defined in two dimensional spaces. In each of the aforementioned 
huhhling disks, a partial hlowup phenomenon occurs, which means if the hlowup 
solutions are scaled according to the maximum of both components in the disk, 
only one component converges to a single equation after taking the limit. The phe¬ 
nomenon of partial hlowup is the major difficulty for understanding the profile of 
huhhling solutions for systems. After identifying huhhling disks we use a Harnack 
inequality (Proposition B helow) proved in |j30)l to describe the behavior of each 
component according to its spherical average around each blowup point. Each 
component is called to have “fast decay” or “slow decay” (see the next section for 
definition) based on its behavior. Roughly speaking, among the two components, 
at least one of them has fast decay and the energy of which is determined. The en¬ 
ergy of the other component can be determined by the Pohozaev identity. Since we 
have to use the Pohozaev identity to determine the energy of one component, this 
approach only works for systems of two equations. Once the energy of at least one 
component is determined in each bubbling disk, we can also do the same for bub¬ 
bling disks in a “group” (see |l30l). A group of bubbling disks looks roughly like 
a single disk after scaling. One major new ingredient in this article is to rule out 
the situation that there are only two bubbling disks in one group. For this we shall 
use Eremenko’s work on surfaces with conical singularities to calculate ^ /iR 2 e“ 
where n is a solution of 



where p\,p 2 are two disjoint points in and 7 i ,72 > — 1 . 

We were not able to prove that all the energy types for G 2 Toda systems are 
multiples of 2 without assumption. The reason that some strange numbers appear 
in Theorem ll.3l is because we don’t have good estimates of the energy of Eiouville 
solution with more than two singular points. 

The organization of this article is as follows. In section two we mention a few 
tools we shall use in the proof of the main theorems. Some majors tools are devel¬ 
oped in OOl and a new tool is based on Eremenko’s work. Then in sections three 
and four we prove the concentration theorems and the proof of a priori estimates 
can be found in section five. 

Acknowledgement Part of the work was finished when fhe second aufhor was 
visiting Taida Insfifufe for Mafhemafical Sciences (TIMS) in December 2014. He 
would like fo express his deep grafifude fo TIMS for fheir warm hospifalify and 
financial aid. 


2. Preliminary discussions 


In fhis section we lisf a number of fools we shall use in fhe proof of Theorem 
11.11 and Theorem 11.31 Many of fhem come from fhe previous work of fhe aufhors 
and J. Wei |[30l and a resulf of A. Eremenko ifT^ . 
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In ll^ Lin-Wei-Zhang studied the concentration of energy for blowup solutions 
to 


(2.1) Aui + '^aijhje^j =47iYi5o, / = 1, ..,n, in B(0,1) 

j 

where A = A„ is the Cartan matrix of order n, are positive smooth functions 

on B(0,1), y; > — 1 indicates the strength of the Dirac mass at 0. Here we recall 
Proposition 2.1 of l(30ll which was the result of a selection process ( see 1381): 

Proposition A: Let be a sequence of solutions to d2. il) with ji = 

■ ■ = Yn = 0. Suppose there is a uniform bound for ijhfe“‘ and the oscillation 

ofu^ on dB{0, 1). Then ifO is the only blowup point ofu^ = {u\,....,u^), there exist 
finite sequences of points := {xj, ....,xf„} (all 0,j = and positive 

numbers l^, .■.,l^„ —)• 0 such that the following four properties hold: 

(1) max;g/{Mf(x*)} = max^^_^ for all j = l,..,m, where I = {l,..,n}. 

(2) exp{\maxiei{u^{x'^f)})l'j -i-oo, j = l,...,m. 

(3) There exists Ci > 0 independent ofk such that 

u'Kx)+ 2\og dist{x,'Lk) < Cl, VxGB(0, 1), i G I = {I, ...,n} 

where dist stands for distance. 

(4) In each B{x^-,l^) let 

(2.2) v\{y) =u\{Eky + x’^j)+2\ogEk, Ek = e^^^'‘, = max max nf. 

Then one of the following two alternatives holds 

(a): The sequence is fully bubbling: along a subsequence (v\,...,v^) con¬ 
verges in C/^^^(]R^) to (vi, ...,v„) which satisfies 

Avi + ^ aijhje'’: =0, i G I. 

j^i 


lim / y aithlle'"' > 4?!, i G I. 

(b):I = JiU J2C ...U J,„U N where J\,J 2 ,...,Jm and N are disjoint sets. 
N and each Jt (t = \ ,..,m) consists of consecutive indices if it has more 
than one index. For each i G N, tends to — 0 ° over any fixed compact 
subset o/M^. The components of v^ = (vj,..., v^) corresponding to each Ji 
(I = I, ...,m) converge in C;^^^(M^) to a SU (|7/| + 1) Toda system, where |7/| 
is the number of indices in 7/. For each i G Jt, we have 


lim 

A:—>00 



> 471 . 


The selection process of Schoen singles out a finite number of bubbling disks for 
the sequence of blowup solutions u^ = {u\,...,u^fj. In each of the bubbling disks, 
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at least one component has energy greater than 2: 



Here we use the integral above to denote the energy of in B{xk,lk)- Since there 
is a uniform bound of the energy of all components, there are only finite bubbling 
disks. The selection process also determines that 


u\{x) + l\ogdist{x,l.k) < Cl, ^ £ fi(0,1) 


which provides a control on the upper bound of the behavior of blowup solutions 
outside the bubbling disks. 

The following proposition in If30ll plays an essential role in the proof of main 
results in their article: 

Proposition B: (Lemma 2.1 of 1(3^ ) For all xq E B(0, 1) \ 'Lk, there exists Co 
independent ofxQ and k such that 

\u^{xi) — uf (a 2 )| < Co, Vxi,X 2 E B{xo,d{xo,Lk)/2), for all i E /. 

Proposition B is a Harnack type estimate which reveals important information 
on the behavior of blowup solutions away from the bubbling area. Let x^ E and 
Tk = jdist{xk,'Lk \ {xk}), then for x,y ^ B{xk, Tk) and |x — = |y — x<;| we have 

u\{x) = Ux^,i{r) + 0{\) where r = |a:<; — 2 c| and 



In other words, the behavior of outside the bubbling disks can be represented by 
its spherical average in a neighborhood of a point in Lk. 

For each Xk GLk, let 



we have 



For each component i we say nf has fast decay on a: E B{0, 1) if 

+2logdist{x,'Lk) < —Nk 


for some Nk If there is a C E M independent of k and 


u^(x)+ 2logdist(x,'Lk) > —C 


we say uf has a slow decay at ;c. Here we note that 

wf (x) + 2logdwt (x,Lk) <C 


holds for all .r E B(0,1). 
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The definition of fast and slow decay is very important for evaluating Pohozaev 
identities. For example in B{xk,r), the following Pohozaev identity holds for solu¬ 
tions to 

” k 

Aui + ^ = 0 

i=\ 

with coefficient matrix being the Cartan matrix (in fact as long as {aij) is symmetric 
and invertible the following holds as well): 


I 

i 

(2.3) = r 






(vV/jf)e“-+2£ 

B{xk,r) I . 

[ +r [ £ {a‘jdyu'^dyu) - \a‘jyu'^yu]) 

JdB(xk^r) ; JdB(xk.r) ; ; ^ 


In order to evaluate the energy concentration from (12.3b it is important to choose r 
so that all components on dB{xk,r) have fast decay. Otherwise the first term on the 
right hand side of (12.3b is not o(l). 

Finally we list a major new tool on the total energy of the following equation: 


r AM + 2e“ =47r(6i — l)5pi +471(02 — l)5p2, in 

(2.4) 

Based on a theorem of Eremenko lIT^ we shall use the following theorem: 

Theorem B-1: Let pi,P 2 be two distinct points in and 0 i ,02 be positive 
integers, then any solution u of M.4^ satisfies 


1 

2n 


e" = 01 + 02 + 03 — 1 


where 03 is a positive integer such that 0 i + 02 + 03 is odd and Oi + 0 y > dk for 
{i,j,k) being any permutation o/(l,2,3). 


Theorem B-1 can be found in the appendix. 


3. Proof of Theorem o 

First we observe that (11.7b is a special case of SU (4) Toda system by letting 
M 3 = Ml in SU (4) Toda system. Consequently Proposition A can be applied to (11.7b 
to obtain the following blowup set = {xj,...,x* } and G M such that all 

x\ (1 = 1, ...,m) tend to the origin and all if (i = I, ...,m) tend to 0. Moreover the 
following properties hold: 

( 1 ) max,{Mf (x^)} = max^g 5 (^ ;p{M 5 (x)}, for all j = 1 , ...,m. 

(2) exp{\m&Xi{u\{x^j)}l]^oo^ j = l,...,m. 

(3) There exists Ci > 0 independent of k such that 

M^(x)+21og(iA?(x,r^) < Cl, VxGB(0, 1), / = 1,2. 

where dist stands for distance. 
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(4) In each B^Xpl'j) let = max,maxg(^^ nf, ek = e 2 ^*^ and 

v’liy) = u-{eky + x^j)+2\ogek, for eky + x^jeB{xpl’^) 

Then one of the following two alternatives holds: 

(a): The sequence is fully bubbling: along a subsequence {v\,v\) con¬ 
verges in to (vi,V 2 ) that satisfies 

r Avi+2e'’‘=0, 

< Av 2 - 2 e^‘+ 2 e '’2 = 0 , in 


lim — 
I :—271 



> 2 , 


/= 1,2 


where (a,y) is the B 2 Cartan matrix. 

(b): = {v\,v\) has only one component converging to a single equa¬ 

tion. If it is the first equation the limit equation is 

Avi-h 2 e'’‘= 0 , in 


If it is the second equation it is 

Av 2 -|- 2 e''^ = 0 , in M^. 


In either case, the convergent component satisfies 

lim -J- [ 2h\e'^> > 2 . 

*- 4=0 271 Jb{x^,i'[) 

Here for convenience we assume lim^,^oo /if (0) = 1. This assumption is not 
essential. 


Remark 3.1. if can be chosen so that for xf G let tt = dist{xp'Lk \ {x^j})> then 
tk/lj °°- 

Since the B 2 Toda system is a special case of the A 3 Toda system, the Pohozaev 
identity for A 3 Toda system can also be applied for the B 2 Toda system. Here we 
recall that for A* = {afj)j,x 3 being the A 3 Cartan matrix we have 

£ a*jOiOj = AY^Oi. 

i.j=i 

Replacing Os by ai we have 

(3.1) 2af-20^02 +02 =401+202- 

Next we consider the energy of global solutions. By the classification theorem 
of Lin-Wei-Ye |[29l . if « = {ui,U 2 -,ui,) is a global solution of SU{3) Toda system 
with finite energy, then 

3 f 

Yaf e+=S7i, I = 1,2,3 
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where {a*j) is the A 3 Cartan matrix. Therefore in our case we let mi = M 3 and then 
we have 

which gives ai = 6 , 02 = 8 . 

It is also easy to verify that ( 6 , 8 ) is the “largest energy” because if (6 + 1,8 + 5 ) 
also satisfies (13.11 ) with s,t >0. Then s = t = 0. 

To understand the concentration of energy we start from any fixed member of 
Say,Xj G Zj. and lef = distance{x\,Lf^\{x\})/2. Let 



Direct computation shows for r G (0, T\^k) 

^ - r \ + (y2,kir) 

^-’ 

dr r 

d-,. 20\,k{r)-202.k{r) 

—U2[r) = -^^-. 

dr r 

The selection process guarantees that 

u\{x)+2\og\x — J{\<C, |x —All < 1 = 1,2. 

If both components have fast decay on dB{x\,r) (r G (0,Ti,^)), {oi^k{r),02,k{r)) 
satisfies 

(3.2) 20i,kir)^ -20i,kir)(y2,kir) + oli^ir) = 40i^k{r) +202,kir) +o{l). 

If we write (13.21) as 

ai,k{r){2ai^k{r) - CJ2,kir)-4) + a 2 ,kir){ 02 ,k{r) - ai^k{r)-2) = o{l), 

we see that for any r, if both components have fast decay, either 2ai ^(r) — 02 ,k > 
4 + o (1), which means 

d . . 2 + 0(1) 

— (Mi(r)+21ogr) <- 

dr r 

or 02 ,k{r) — 01,(1 > 2 + 0(1), which implies 

d . _ . \ . 2 + 0(1) 

— (M 2 (r)+ 21 ogr) <-. 

dr r 

If one component (say u\) satisfies 

^(Mi(r)+21ogr) >0, 
dr 

there is a possibility that for some larger radius 5 , u\ becomes a slow decay com¬ 
ponent on dB{x\,s). 

Next we consider the possible energy concentration types in B{x\,'Zi^k)- Let 
B{xf[,l\) be a bubbling disk. Suppose in this disk, the first component converges 


f 201-02 = 4 , 

\ - 201+202 = 4 
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according to the scaling of the maximum. If the region is close enough to the 
center of blowup ( the region still tends to infinity after scaling according to the 
maximum), we have 

0 'Lt:( 4 ) = 2 + 0(1), 02^k{h) = o{\). 

Now we consider the energy change from B{x\,lk) to First we notice 

that on dB{x\,lk), ^(^2 + 21ogr) > 0, which means u\ may become a slow decay 
component when r increases. So the first possibility in B{xf[,Ti^k) is that 

aLI:(Tl,I:) =2 + 0 ( 1 ), O2,k{'lhk)=o{l), 

which means u\ does not change to be a slow decay component. It is proved in 
OOl that if no component changes to a slow decay component, the energy of each 
component only changes by o (1) ( see Lemma 5.1 of . Even though that lemma 
addresses SU (3) Toda system but a very similar proof also applies to this case). 

Since ^(w 2 + 21ogr) > 0, u\ could become a slow decay component before r 
reaches Xi^k- Suppose at some s> r, 

u\{s) +21og5 > —C 

for some C > 0 very large. Here we observe that at this moment starts to in¬ 
crease its energy but the energy of u\ barely changes because + l log r) is 

still negative. which means Ti ^ is very large comparing with s, there 

is N such that at dB{ 2 J\,Ns) 

02 ,k{Ns)>5, Oi^k{Ns)=l + o{\) 
u\{Ns) +l\og{Ns) < —Nk, for some Nk^°° 

^{u2 + 2logr)\r^N, < 0, ^(mi +21ogr)|^=A,^ > 0. 

In other words, from r = s to r = Ns, the energy of increases and as a result, 
the derivative of U 2 {r) + 21ogr changes from positive to negative. But because of 
the Harnack inequality (Proposition B) Mj is still a fast decay component and its 
energy barely changed, even though at r = the derivative of ui +21ogr has 
become positive due to the change of the energy of u^- Since Ti^k/^ —)• 00 we can 
find fending fo 0 slowly such fhaf Nj^s < X\^k/2 and on dB{x\,Nls) bofh u\ and 
have fasl decay. Evaluafing fhe Pohozaev idenfify on dB{x^^,N'j^s) we have 

ai,,t(iV(^)=2 + o(l), 02 ,kiN'ks) = 6 + o{\). 

If Ti^k is only comparable fo s, fhen on dB{x\,T\^k)^ “2 ^ ^1°"' decay componenf 

and Oi^k{xi,k) = 2 + o{\). 

At dB(x\,N',^s), 

^(u1(r)+2logr) = i±^, r = N',s, 
dr r 

and 

^(M|(r)+2logr) = — r=N'j,s. 
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So at this radius there is a possibility that u\ may become a slow decay component 
for larger r. 

By exactly the same reason as before it is possible that u\ increases to a slow 
decay component on for which the second component has the energy 

=6 + 0 ( 1 ); or Mj finishes its transition of energy before r reaching X\y. 
such that both components have fast decay on and 

(y\,k{sk) = 6 + 0 ( 1 ), 02 ,k{sk) = 6 + 0 ( 1 ). 

Similarly if Sk is small compared to we could also have ( 6 , 8 ) as the energy 
type in B(xj, Ti,jt). Since ( 6 , 8 ) is the type of the energy a global solution has, there 
is no extra energy outside (see Theorem 4.1 and Theorem 4.2 of 1301 ). 

If we start with the type ( 0 , 2 ), which means in B{x\,lk), <^ 1 ,^( 4 ) = o(l) and 
<y 2 ,k{h) =2 + 0 ( 1 ). Then the following type may occur: 

(3.3) (0,2), (4,2), (4,8), ( 6 , 8 ). 

More specifically if bofh componenfs have fasl decay on dB{x\,Ti^k) then 

{(yx^k{'Ci,k),02,k['^i,k)) = (a + o(l),Zi + o(l)) 

where {a,b) is one of fhe four types in (13.31) . 

If one of fhe fwo componenfs has slow decay on dB{x\,T\^k)^ from fhe dis¬ 
cussion above we see fhaf fhe energy of fhe ofher componenf (which has fasl de¬ 
cay) is a mulliple of 2. For example, if fhe second componenf has slow decay on 
5 B(xj,Ti^<:), is 2 + 0 ( 1 ), 4 + 0 ( 1 ) or 6 + o(l). If fhe firsl componenf has 

slow decay, 02 ,k{'^\,k) is 2 + o(l), 6 + o(l) or 8 + o(l). 

Now we consider bubbling disks in a group. The concepf of group is infro- 
duced in l30l ). which means we consider bubbling disks relatively close fo one 
anofher buf relatively far away from ofher members in For example, if Xj 
are in one group, dist{x^,x'j) ~ dist{x^,x^) for i,j,m,n G {1,2,3} buf i / j and 
n (here means comparable). Moreover for any G buf a 0 {1,2,3}, 
dist{x^, 2 J[)/dist{x^,x^^) —)• 00 . By Proposition B if bofh componenfs have fasl de¬ 
cay around one bubbling disk, bofh componenfs have fasl decay around any of fhe 
disks in Ihis group. Suppose fhe group members are B(x}, ty^^,), ...,B(x^, Tm,k)- By 
fhe definilion of group, all fhe are comparable. If bofh componenfs have fasl 
decay we can find Nk^°° slowly such lhal all members in Ihis group are conlained 
in B{x\,NkXi,k) and 

m 

(yi.ki4^^kT:i,k) = £ (yi,ki^pT:j,k)+o{l), 

1=1 

m 

02,k{A^NkT:i^k) = Yj t^2.;t(Y/,Ty.,t)+o(l), 

1=1 


where 


0\,k{x\,Nk'Ci,k) = 


2n 


h\d 


and ofher nolalions are underslood similarly. 
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In Other words, Nk can be chosen in a way that both components still have fast 
decay on dB{x\,Nk'C\^k) and the energy of each component in this larger region is 
a o(l) perturbation of the sum of the energy in each bubbling disk of this group. 
As observed before around each bubbling disk, at least one component has fast 
decay and the corresponding energy is a multiple of 2 + o(l). Therefore both 
Oi.k{^\.,Nk‘X\^k) and 02 ,k{^\,NkTi^k) are multiples of 2 and they must satisfy the Po- 
hozaev identity. The following are the only pairs that satisfy the Pohozaev identity 
with each component being a multiple of 2: 

(2,0), (0,2), (4,2), (2,6), (4,8), (6,6), (6,8). 

The final case we consider is when only one component has fast decay in a 
group. Suppose u\ has fast decay in the group described as before. In this case 
is a multiple of 2. In each bubbling disk, say '^j,k), since 
has slow decay, 02,k{xp '^j,k) > 0 and by Proposition B we can choose Nk^ °° so 
that u\ is still a fast decaying component on dB{2i[,Nk'Ci^k) and 

m 

Oi.kix\,Nki:i^k) = £ Oi^kix'j,T:j^k)+o{l). 
i=i 

Moreover, also has fast decay on dB{x\,NkTi^k)- Therefore the Pohozaev iden¬ 
tity can be evaluated on this radius. Obviously this group is contained inB(vj ,A<.Ti ^) 
Since the first component is a multiple of 2, we see immediately that there is no 
new type except those we have known. For example, if we have two regions that 
both grow up from the type (2,0), then the energy of the second component on 
B(jc^,A<:Ti^<:) has to be 8 + o(l). If the group has two regions that grow from (2,0) 
and (4,2), the energy of the second component may grow to 8 and make the en¬ 
ergy type (6,8). If there are types of (2,0) in this group, the energy of on 
B{x\^NkT\^k) has to be 6 -|-o(l) to 8 -|-o(l). If the energy of u\ in this group is 
greater than 6, for example 8 -|-o(l), we use scaling to make the distance between 
any two blowup points in this group comparable to 1. Then u!^, since it has slow 
decay in the neighborhood of these points, converges to a function U 2 that satisfies 

m 

(3.4) Au 2 + 2e‘‘^ =47i{Y^Yj5pj), in 

./=i 

where pj {j = are the limits of blowup points in the group after scaling. 

Each Yj is a multiple of 2 and Yj = By standard potential analysis it is easy to 
see that there exist a > 2, C G M and a > 0 such that 


M2(-^) = —Ctlog|x|-|-C-|-o(|x| ^), |x| > 1 


and 

Vu2{x) = — a-j-^ -h 0{\x\^^^^). 
Integrating both sides of (13.41 ) we see that 


1 

2n 





>9. 


(3.5) 
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This estimate means in the group contained in if = 

8 + o(l), 

02 ,k{AMT:u) >9. 

It is easy to see from direct computation that this pair of numbers cannot satisfy 
the Pohozaev identity. 

If has fast decay in the group mentioned before, and u\ has a slow decay 
in the neighborhood of the aforementioned bubbling disks, we first remark that if 
the total energy of u\ in these bubbling disks is 6 + o(l), the energy of u\ in, say, 
is either 2 + o(l) or6 + o(l), according to the Pohozaev identity. If 
the total energy of u\ in this group is 8 + o(l), the total energy of u\ in 
is either 4 + o(l) or 6 + o(l). If the total energy of u\ is 2m + o(l) for m > 5,m G N 
(the set of natural numbers), as in the previous case we first derive a lower bound of 
the energy of u\ in By scaling the distance between any two members 

of the group is comparable to 1, then the slow decaying component u\ converges 
to Ml, which satisfies 

(3.6) Ami+ 2e"'= 27r2^7,-5p., in 

j 

where pj are the images of blowup points in after scaling, 7 y = 2m for m>5. 
By standard potential analysis 

Vmi(. r) = —a-j—p + l-icl > 1 

for some a >2 and a > 0. Integrating both sides of (13.61) we have 

— [ > m + 1. 

2k Jr2 

Direct computation shows that there exists no pair (ai, 02 ) satisfying the Pohozaev 
identity with 02 = 2m (m > 5) and ai > ni + 1. 

Finally we rule out the case that there are two (0,2) type bubbling disks in the 
group and they are the only members. Note that this is the only case that the energy 
of m| can be 4 + 0 ( 1 ). Since u\ is slow decay. We see that by scaling the distance 
between these two groups into 1, we see that the re-scaled version of u\, which we 
use M to represent, satisfies 

Am + 2o“ = 471+ 4715p2 in 

and we also know fj ^2 0 " < o®. By Theorem B-1 we know the total integration of 

J— [ 0 " = 4 or 6. 

2k Jr2 

But neither (6,4) nor (4,4) satisfies the Pohozaev identity, which means it is not 
possible to have two (0,2) type regions in one group. 

Since the combination of groups is similar to those of bubbling disks in one 
group, we have exhausted all the concentration types. Theorem 1 1.1 1 is established. 
□ 
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4. Partial results for G 2 Toda system 
The G 2 Toda system we consider is 

Aui + 2/iie“* — /i2e“^ = 0, 


'• AMi-3/iie“>+2/i2e“2 =0, in B(0,l)c: 

and (11.8b and (11.9b also hold for as a sequence of solutions to (14.1b . 

The G 2 Toda system is a special case of a Ae Toda system where = U 4 = 
u\ +log2, M 5 = U 2 and U(, = ui. Using the Pohozaev identity for Ag one obtains 
easily the following identity: 


(4.2) 


3 ai^ + a| - 3 ai 02 = 6 ai + 202 - 


to 


First we consider global solution. We apply the theorem of Lin-Wei-Ye 
have 

6 

e‘'J = Sn, / = 6 

where (a,y)6x6 = Ag. Using 03 = 04 = 2ai, 05 = O 2 , <^6 = we see the blowup 
type for global solution is (12,20). 

Next we consider partial blowup cases. Starting from B{x\,T\^k) which comes 
from the selection process as before. First we assume 

O'!,it( 4 ) = 2 + 0(1), 02,k{h) = o{\). 

From here we consider all the possible concentration types as r increases from 4 
to Ti^k- If k -2 becomes a slow decaying component the energy type could change to 
(2,8). After this (8,8) could occur and finally (8,18) could occur. If we start from 
(0,2), then u\ may change to a slow decaying component that leads to (4,2). From 
here we could have (4,12) and (10,12). 

So far the following types are possibilities: 


(4.3) (2,0),(0,2),(2,8),(4,2),(12,20), 

(4,12),(8,8),(8,18),(10,12). 

When we consider the derivative of spherical averages, we have 

d_ -20u{r) + 02,k{r) d_ 3ai^k{r)-2a2,k{r) 

—ui[r)- -, —U2(r)- -. 

dr r dr r 

If Pohozaev identity can be evaluated on B{x,r), which means both u\ and are 

4 + 0 ( 1 ) 


4 + 0 ( 1 ) 


fast decaying 

on dB{x, r), either we have 


(4.4) 

, , 4 + 0 ( 1 ) 

u\{r) < -, and 

r 

U2{r) > 

or 



(4.5) 

- ^ N, ^ 4 + 0 ( 1 ) 

ui[r)> -, and 

U2{r) < 


On the other hand if both 








TODA SYSTEM 


17 


(4.6) loi^kir) - <y 2 ,k{r) > 2, and -3ai^k{r)+2a2,kir) > 2 

on a fast decaying radius for both components, then there is no essential energy 
outside. 

Next we consider the combination of bubbling disks in a group. If both compo¬ 
nents have fast decay in the neighborhood of those bubbling disks, both 0 \^k{x\,NkT\^k) 
and (y 2 ,k{^i,Nk'^\,k) are multiples of 2 and all such pairs that satisfy the Pohozaev 
identity have been listed in (14.31) . So we only consider the case that one component 
has slow decay in the neighborhood of those bubbling disks in one group. If the 
group has only two disks, we shall use Theorem B-1. From Theorem B-1 we see 
that the energy of the slow decaying component is also a multiple of 2. So we still 
do not add any new type except those listed in (14.31) . 

Finally we rule out the case that there are at least three bubbling disks in one 
group. The most “energy efficiency” case is there are three (2,0) type bubbling 
disks in one group or three (0,2) type bubbling disks in one group. In the first 
case, clearly has slow decay and it is easy to see that 

(6,10 + 2^7) 

satisfies the Pohozaev identity. By assuming 02 < 10 + 2y/l we ruled out this 
case. If u\ has the slow decay and there are three (0,2) type bubbling disks in the 
group, we see easily that (4 + 2\/2,6) satisfies the Pohozaev identity. By assume 
01 < 4+2\/2 this case is also ruled out. Other cases of having more bubbling disks 
or only three bubblings disks with more energy can also be ruled out easily by the 
restriction of 0,. 

Theorem [T3]is established. □ 

5. Proof of Theorem 1 1.2 1 and Theorem 1 1.41 

Since the nature of these two theorems is so close we just prove Theorem [T2] as 
an example. The proof of Theorem 1 1.4h s completely similar. 


Let 



Ui = Ui — log / , / = 1,2. 


Jm 

Then we have 

( AgMi +2pi(/iie"‘ - 1 )-P 2 {h 2 e"^ - 1) =0, 

(5.1) 

1 AgU 2 — 2pi(/iie"‘ — 1) + 2p2(/i2e"^ — 1) = 0, 

and 


(5.2) 

f hie"'dVg = l, / = 1,2. 

Jm 


In the first step, which is the major one we prove that there is a C independent 
of u such that 

Ml I + |m2| < C. 


(5.3) 
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In order to prove (15.31) we only need to prove 


(5.4) max max M,- < C. 

i M 

because once (15.41) is established the lower bound can be obtained by standard 
Hamack inequality. Thus we shall prove (15.41) by way of contradiction: Suppose 
= {u\,U 2 ) is a sequence such that maxjmax^u^ —)• oo. Let G be the blowup set 
for u^: 

G = {p€M; —)• p, limmaxu^{xk) ^ } 

>oo i 

Let pi,...,pN be the blowup points of = {u\,u^ onM. It is easy to prove that 
there are only finite blowup points on M by standard estimates ( see |l7j|33]|28l, etc). 
Here we assume that pi,..,pyv are distinct points. By the Green’s representation 
of u'l it is easy to see that u'l has bounded oscillation outside the bubbling area: 
U^iB(p,', 5()) where 5o > 0 is chosen small enough to make 

B(p,-,5o)nB(py,5o) =0. 


Around each p/, we write the equation in Euclidean form in local coordinates, 


let 


1 


a,(pi) := lim lim — 


plhie“-dVg, 


5 —> 0 ^:^“ 2,71 Jb{p[,S) 

then we shall see that this is exactly the same as what is defined for fhe local 
equafion. Indeed, lef p\ —>• p\ be where max,maxg(pj iif is attained. Take fhe 

L 7 ¥{y k) . -J Ox, 

local coordinates around p\, fhen ds has fhe form e ''i [dyf + dy^) where 
I Vva( 0)| = i/r(0) = 0, Ay= -IKe'f', 

K is fhe Gauss curvafure. Obviously fhe equafion for can be written as 
+ 2pfe’^/iie"^ +e'^(-2pf+ p|) = 0, 


(5.5) 


Aii^ - 2pfe'^/iie"i + 2p^e'i'h2e"2 +eV'(2pf - 2p|) = 0, 

Lef /f and /| be defined by 

A/f + r.'^(-2pf + p|)=0, in Bg,, 
on dBs^. 


and 


f =0, 


A/| + r.v^(2pf-2p|)=0, 

on dB^ 


in B 


SqI 


/f= 0 , on ODSg. 

By seffing = pfhie'^^f> and u\ = we can write (15.5b as 

Au\ + 2h\e^' - 11 ^ 6^2 = 0, in 


in B 


5o 


Au\-2h\e‘^i +2/i|e“2 =0, 

Since dVg = e'^dy if is easy fo see fhaf 


h^e^'dy = 


Ib(p\,So) 


pfhie°idVg, 


i= 1,2. 
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Then (<7i(pi),0'2(pi)) is just one of the seven types. Moreover, at least one 
of the components has fast decay. For example if (ai(pi),a 2 (pi)) = (4,2), the 
second component decays fast. Since the oscillation of either component is finite 
away from the bubbling disks, we see that at least one component has little energy 
outside the bubbling disks. Also if the second component is fast decaying, the 
second component around any blowup point is also fast decaying. Then we have 


U 


N 

t=\ 


2n 



for some positive integer N. However since 


Im\ub(p,,S) 


p^h2e"^dVg = o(l) 


and (15.21) holds, we get a contradiction to our assumption that p| cannot tend to a 
multiple of An. Thus (15.41) . and consequently (15.31) are established. 

Now we finish fhe proof of Theorem 1 1.21 Clearly we can wrife fhe equafion for 
u = (mi,M 2 ) as 

{ higUi +2pi(/iie“‘ — 1) — P2(/i2«“^ — 1) = 0) 

AgU2 — 2pi (h\e^^ — 1) T 2p2ih2B^^ — 1) = 0. 

Since 


[ f2pi(/qe“‘ 
Jm V 

[ f-2pi(/iie“‘ 
Jm V 


-\)-p2{h2e^^-\)^dVg=0, 
-l) + 2p2ih2e^^-l)]dVg = 0 


and (15.31) holds, from sfandard elliptic esfimafe and = 0 (/ = 1,2) we see fhaf 
|ui I + |m 2 | < C. Theorem 1 1.21 is esfablished. The proof of Theorem 1 1.41 is similar. 
□ 


6. Appendix: Two Theorems of Eremenko 

In (his appendix we inferpref some fheorems of A. Eremenko |[T9ll info a form 
thaf can be used in fhis arficle. 

Eirsf we recall fhaf on a Riemann surface S, a mefric go is called conformal if in 
any local coordinafe sysfem z/ G C C, 

goizi) =e‘'\dzi\^, zi^O. 

for a measurable and bounded function u in Q.. Eel S, 2 i conformal mefric gs 
is called (o have a conical singularify a( Pa of (ofal angle 2n{a + !)(«> —1) if 
(here exisf local coordinafes z{P) G C C and u G C*’(fl)nC^(fl\{P}) such fhaf 
z(Po) = 0 and 

where gs is fhe local expression of g^. 
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In lfT9]| A. Eremenko studied the following situations: Let pi,p2,P3 be distinct 
points on and iTidjij = 1,2,3) be their total angles of conic singularity respec¬ 
tively {Gi > 0 ,i= 1,2,3). 

The first case is at least one of the three numbers G\ ,62,63 is an integer. Among 
other things Eremenko proved: 

Theorem A (Eremenko): Let 6\> \ be an integer, but 62 > 0 and 63 > 0 are not 
integers. If there is a conformal metric with curvature 1 on the sphere, with three 
conic singularities of angles ItlGi (i = 1 , 2 , 3 ), then either 62 + 63 or {62 — 63\ is an 
integer m of opposite parity from 61, and m<6\ — \. 

Remark 6.1. Two integers A and B are called to belong to opposite parity if one 
of them is even and the other is odd. 

Theorem A is the existence part of Eremenko’s original statement in |[T9ll . the 
uniqueness part of conformal metric with prescribed conical singularities can be 
found in Eujimori, et. al 1201 . 

If all 6i are positive integers not equal to 1, Eremenko proved: 

Theorem B (Eremenko): Let 61,62,63 be three positive integers not equal to 1. 
If there exists a conformal metric of curvature 1 on the sphere, with three conic 
singularities of angles 1 ti6\, 27102 tmd respectively, then 0+ 02 + 03 is odd 
and 0, + dj > 6k for {i,j,k) being any permutation of ( 1 , 2 , 3 ). 

The second case is none of 61,62, Gt, is an integer. This is case is not used in this 
article but we still translate it into a PDE result for applications in the future. 

We say (0i ,62,63) ( each 0, > 0) is equivalent to (±0i + m, ±02 ± n, ±03 ± k) 
when (m,n,k) are integers with the property m + n + k = 6 ( mod 2). Every non¬ 
integer triple is equivalent to one and only one triple with the property 

(6.1) O<0i'±02<l, O< 02 ± 03 < 1 , O<0(±03<1. 

Eor this case Eremenko proved: 

Theorem C (Eremenko): If none of 61,62,63 is an integer, then a conformal 
metric of constant positive curvature on the sphere with conic singularities of total 
angles 2 ti 6 i, 27162 and 27163 exists if and only if the unique equivalent triple with 
the property d6. iP satisfies 6[ + 62 + 6^ > 1. Such a metric of curvature 1 is unique. 

We shall interpret Theorems A,B and C for the following equation: 

(6.2) An ±2e" = 471(01 -l)5„i± 471(02 -1)5„2, in f e " < 00 . 

where ,p 2 are two distinct points in and 0i , 62 are positive constants. 

Theorem A- 1 : Let pi, p2 be two distinct points in R^ and ube a solution o/Q. 
Suppose 01 is a positive integer and 62 > 0 is not an integer, then any solution u of 
/ 16 . 2 I) satisfies 


± [ e" = 0i±02±03-l 
271 Jr2 
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for some 63 > 0. Moreover either 62 + 63 or 1 62 — 63 ! an integer m of opposite 
parity from d[, and m< 61 — 1 . 

The following theorem is a translation of Theorem B: 

Theorem B-1: Let pi,P 2 be two distinct points in and u be a solution of 
f l 6 . 2 l) . If both 61 and 62 are positive integers, u satisfies 

± f e" = 61 + 62 + 63-1 
271 Jr 2 

where 63 is a positive integer, 61 + 62 + 63 is odd and 6 i + 6 j > 6 k for (/, j, k) being 
any permutation o/(l,2,3). 

The interpretation of Theorem C (which we don’t use in this article) is 

Theorem C-1: Suppose 61 , 62,63 are all positive but none of them is an integer. 
Then there exists a unique solution u of ^6.2^ with 

— [ 0 “ = 01 + 02 + 63 — 1 

271 Jr 2 

if and only if 6 [ + 02 + 03 > 1. 

Proof of Theorem A-1,Theorem B-1 and Theorem C-1: 

Step one 

In the first part of the proof, we state one fact: Let qi, q 2 and ^3 be distinct points 
on C and 2710,-, (/ = 1,2,3) be the total angle of conical singularity at qi. Then the 
equation for Gauss curvature equal to 1 is 

(6.3) Av + e^'’ = 27i{6i - 1)5^, +27r(02 - 1)5^^+27r(03 - 1 ) 5 ^ 3 . 

Here is the reason why (16.31) holds. First away from the singularities 

K = - {Av)e-^'’ = \ . 

Then around each singularity, say qi, if in the neighborhood we let 

V = V - (01 - l)log|Y-^i|, 
the equation for v around <71 would be 

+ =0. 

Thus 27101 is the total angle at qi. Let vi = 2v, then (16.31 ) becomes 

(6.4) Avi +2e^’' =47r(0i - 1)5^, +471(02-1)5^3 +47r(03 - 1)5^3. 

Theorem A and Theorem B can be applied if 00 is a not a singular point. 00 is 
not a singular point if 

vi(y) = -41og|y|+0(l), |y|>l 

because if we let V 2 (z) = vi(z/|z|^) — 41og|z| for z close to the origin, we would 
have 

Av2(z) + = 0 

in a neighborhood of 0 and V2 is bounded near 0. 
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Step two: we can assume p\ = (—1,0), p 2 = (1,0): 

Going back to equation (16.21) . we first use a transformation to move p\,p 2 in 
(16.21) to (—1,0) and (1,0): Let 

ui{x) = u{x)-l{Bi - l)log|v-pi| - 2(^2 - l)log|x-;i 2 |, 
then clearly 

Ami(x) = Am- 471(01 - l)5p,- 471(62 - l)5p2, in M^. 

and 

Ami+2c“'|v-;7i|2(®‘-i)|;c-P 2|^^®-^^^ =0, in 

Let 

uiiy) = ui{ ^^ +r/y)+2(0i + 02 - l)logr/, d = \p\-p 2 \- 

Then it is easy to verify 

AM 2 (y) + 2 |y-Ap^^‘^^^|J-^ 2 p^®"^^^^’"" =0, in 
where Pi = (—1,0), P 2 = (1,0). Clearly 0 is a regular point for U 2 - 

Step Three: Singularity at infinity 

Let 

U 2 {y) =M2(y)+2(0i - l)log|y-Pi| +2(02- l)log|y-P 2 |, 

we have 

Am3 = Au2 + 47r(0i — 1) + 471 (02 — 1) 5p2 

and 


(6.5) Am3+2c" 3 =4;r(ei - l)5p^ +471(02-1)5^2. 

Then we consider the Kelvin transformation of M3: 

M 4 (z) = M 3 (- 4 r-) -41og|z|, zGM^, 

zr 


then round 0, 


Am4(z) +2e“‘* = 0, 


m 


fil/2\{0}. 


In order to determine AM4 (0) we first observe that by standard potential analysis 
M3(y) = -alog|y|+ci + 0(|y|^^), |y| > 1, 


(6.6) 


VM3(y) = -a^ + 0(|y|—1) 


bl>i- 


' |y| 

for some a > 2, ci G M and a > 0. 

By the definition of M3 and M4 we see that (16.61) leads to 

M4(z) = (a — 4)log |z| + 0(1) near 0 


and 

oc 

Am4(0) = 2n{a — 4)5o = 47r(y — 2)5o. 


If we use 27103 to denote the total angle at 0 we have 
(6.7) = 
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Step four: Completion of the proof: 

Now we integrate on both sides of (16.51) . the left hand side gives (using (16.61) ') 

lim / Am 3 = —Ina = —271(203 +2) 

R^°°Jbr 

where the last step is by (16.71) . Direct computation shows that 



Thus 


-271(203+2)+2 [ 471(01-1)+471(02-1), 

and we have obtained the following important equation: 


( 6 . 8 ) 


-f 

271 Jm? 


d[ + 62 + 6^ — 1. 


To finish the proof of Theorem A-1, we see that it follows directly from Theorem 
A if 01 / 1. If 01 = 1, it is easy to verify that Theorem A-1 still holds, because the 
requirements on 62 and 03 imply 02 = 03 and it is easy to check that the conclusion 
still holds by Prajapat and Tarentello’s classification theorem llTTl . See 1401 as 
well. Theorem B-1 follows from Theorem B if no 0, is 1. If some 0, is one, say 
01 = 1, the classical result of Troyanov BOl asserts 62 = 63 . It is easy to see that 
the conclusion of Theorem B-1 still holds in this case. Finally for this case we note 
that 01 = 02 if 03 = 1. Theorem C-1 is straight forward interpretation of Theorem 
C. □ 


References 

[1] Ao, W., Lin, C. S., Wei, J. On Non-topological Solutions of the A 2 and B 2 Chem-Simons 
System, Memoirs of American Mathematical Society, in press. 

[2] Bartolucci, D.; Tarantello, G. The Liouville equation with singular data: a concentration- 
compactness principle via a local representation formula. J. Differential Equations 185 (2002), 
no. 1, 161180. 

[3] Battaglia, L; Malchiodi, A; A Moser-Trudinger Inequality for the singular Toda system, 
preprint. 

[4] W. H. Bennet, Magnetically self-focusing streams, Phys. Rev. 45 (1934), 890-897. 

[5] Bolton, J., Woodward, L.M.: Some geometrical aspects of the 2-dimensional Toda equations. 
In: Geometry, Topology and Physics, Campinas, 1996, pp. 69-81. de Gruyter, Berlin (1997). 

[6] Bolton, J., Jensen, G.R., Rigoli, M.,Woodward, L.M.: On conformal minimal immersions of 
S2 into CPn. Math. Ann. 279(4), 599-620 (1988). 

[7] H. Brezis, F. Merle, Uniform estimates and blow-up behavior for solutions of —Am = V (jc)e“ in 
two dimensions. Comm. Partial Differential Equations 16 (1991), no 8-9, 1223-1253. 

[8] Calabi, E.: Isometric imbedding of complex manifolds. Ann. Math. 58(2), 1-23 (1953). 

[9] Ching-li Chai, Chang-shou Lin, Chin-lung Wang, Mean field equations, hyperelliptic curves, 
and modular forms: I, preprint 2014. 

[10] C. C. Chen, C. S. Lin, Sharp estimates for solutions of multi-bubbles in compact Riemann 
surfaces. Comm. Pure Appl. Math. 55 (2002), no. 6, 728-771. 

[11] C. C. Chen, C. S. Lin, Topological degree for a mean field equation on Riemann surfaces. 

[12] Zhijie Chen, Ting-jung Kuo, Chang-shou Lin, Mean field equation, isomonodromic deforma¬ 
tion and Painleve VI equation: Part I, preprint, 2015. Comm. Pure Appl. Math. 56 (2003), no. 
12, 1667-1727. 


24 


CHANG-SHOU LIN AND LEI ZHANG 


[13] Chern, S.S., Wolfson, J.G.: Harmonic maps of the two-sphere into a complex Grassmann man¬ 
ifold. W.Ann. Math. 125(2), 301-335 (1987). 

[14] Doliwa, A.: Holomorphic curves and Toda systems. Lett. Math. Phys. 39(1), 21-32 (1997). 

[15] G. Dunne, R. Jackie, S.Y.Pi, C. Trugenberger, Self-dual Chem-Simons solitons and two dimen¬ 
sional nonlinear equations, Phys. Rev. D 43 (1991), 1332-1345. 

[16] G. Dunne, Mass degeneracies in self-dual models, Phys. Lett. B 345 (1995), 452-457. 

[17] G. Dunne, Self-dual Chem-Simons theories. Lecture Notes in Physics, Springer, Berline, 1995. 

[18] G. Dunne, Vacuum mass spectra for SU {N) self-dual Chern-Simons-Higgs, Nucl. Phys. B 433 
(1995), 333-348. 

[19] A. Eremenko, Metrics of positive curvature with conic singularities on the sphere. Proc. Amer. 
Math. Soc. 132 (2004), no. 11, 3349-3355. 

[20] S. Fujimori, Y. Kawakami, M. Kokubu, W. Rossman, M. Umehara, K. Yamada, CMC-1 trinoids 
in hyperbolic 3-sphere and metrics of constant curvature one with conical singularities on the 
2—sphere. Proc. Japan Acad. Ser. A Math. Sci. 87 (2011), no. 8, 144-149. 

[21] N. Ganoulis, P. Goddard, D. Olive, : Self-dual monopoles and Toda molecules. Nucl. Phys. B 
205, 601-636 (1982). 

[22] Guest, M.A.: Harmonic Maps, Loop Groups, and Integrable Systems. London Mathematical 
Society Student Texts, vol. 38. Cambridge University Press, Cambridge (1997). 

[23] J. Jost, C. S. Lin and G. F. Wang, Analytic aspects of the Toda system II: bubbling behavior and 
existence of solutions. Comm. Pure Appl. Math. 59 (2006), no. 4, 526-558. 

[24] J. Jost, G. F. Wang, Classification of solutions of a Toda system in R2. Int. Math. Res. Not. 
2002, no. 6, 277-290. 

[25] H. Kao, K. Lee, Selfdual S(/(3) Chern-Simons-Higgs systems, Phys. Rev. D 50 (1994) 6626- 
6635. 

[26] K. Lee, Relativistic nonabelian Chem-Simons systems, Phys. Lett. B 225(1991), 381-384. 

[27] K. Lee, Selfdual nonabelian Chem-Simons solitons Phys. Lett., 66 (1991), 552-555. 

[28] C. S. Lin, J. C. Wei, W. Yang, Degree counting and shadow system for SU{3) Toda system: 
One bubbling, preprint 2014. 

[29] C.S. Lin, J. C. Wei, D. Ye, Classifcation and nondegeneracy of SU (m+ 1) Toda system. Invent. 
Math. 190(2012), no.l, 169-207. 

[30] C. S. Lin, J. C. Wei, Lei Zhang, Classification of blowup limits for SU(3) singular Toda systems. 
Analysis and PDE, in press. 

[31] C. S. Lin, Lei Zhang, Profile of bubbling solutions to a Liouville system. Ann. Inst. H. Poincar 
Anal. Non Linaire 27 (2010), no. 1, 117-143. 

[32] C. S. Lin, Lei Zhang, A topological degree counting for some Liouville systems of mean field 
equations. Comm. Pure Appl. Math, volume 64, Issue 4, pages 556-590, April 2011. 

[33] C. S. Lin, Lei Zhang, On Liouville systems at critical parameters. Part I: One bubble. J. Funct. 
Anal. 264 (2013), no 11, 2584-2636. 

[34] A. Malchiodi and D. Ruiz, a variational analysis of the Toda system on compact surfaces. 
Comm. Pure Appl. Math. 66 (2013), no. 3, 332-371. 

[35] Mansfield, P. Solutions of Toda systems. Nucl. Phys. B 208, 277-300 (1982). 

[36] Nolasco, M., Tarantello, G.: Vortex condensates for the SU(3) Chem-Simons theory. Commun. 
Math. Phys. 213(3), 599-639 (2000). 

[37] J. Prajapat, G. Tarantello, On a class of elliptic problems in R^: symmetry and uniqueness 
results. Proc. Roy. Soc. Edinburgh Sect. A 131 (2001), no. 4, 967-985. 

[38] R. Schoen, Stanford lecture notes. 

[39] M. Troyanov, Prescribing curvature on compact surfaces with conical singularities, Trans. AMS, 
324 (1991) 793-821. 

[40] M. Troyanov, Metrics of constant curvature on a sphere with two conical singularities, Lect. 
Notes Math., 1410, Springer, NY, 1989, 296-308. 

[41] K. Uhlenbeck, Harmonic maps into Lie groups: classical solutions of the chiral model. J. Dif¬ 
ferential Geom. 30 (1989), no. 1, 150. 



TODA SYSTEM 


25 


[42] Y. Yang, The relativistic non-abelian Chern-Simons equation. Commun. Math. Phys. 186 (1), 
199-218 (1999). 

[43] Y. Yang, Solitons in Field Theory and Nonlinear Analysis. Springer Monographs in Mathemat¬ 
ics. Springer, New York (2001) 

Taida Institute for Mathematical Sciences, Center for the Advanced Study in 
Theoretical Sciences, National Taiwan University, Taipei 106, Taiwan 
E-mail address: cslinOmath. ntu. edu. tw 

Department of Mathematics, University of Florida, 358 Little Hall P.O.Box 
118105, Gainesville FL 32611-8105 
E-mail address: leizhangSuf 1. edu 



